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DGLAP evolution extends the triple pole pomeron fit 

G. Soye^ 

Inst, de Physique, Bat. B5, Universite de Liege, Sart-Tilman, B 4 OOO Liege, Belgium 

We show that the triple pole pomeron model P| provides an initial condition for a DGLAP 
evolution @ that produces a fit to high experimental DIS data. We obtain good for initial 
scales down to 3 GeV^. Values of the initial scale smaller than 1.45 GeV^ are ruled out at the 90% 
confidence level. 
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I. INTRODUCTION 


We have shown in a previous paper Q that it is possi¬ 
ble to fit the experimental data for F 2 with a double or 
triple pole pomeron model in the region 


2v > 49 GeVG 

cos(0t) = ^ > 


< 150GeV^ 

X < 0.3. 

V — 


49 

2m^ ’ 
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We have also shown that one can extend the usual t- 
channel unitarity relations to the case of multiple 
thresholds and multiple poles. This allowed us to pre¬ 
dict from F^ and the pp total cross-section. In the 
latter case, we have shown that all processes have the 
same singularity structure. 

However, in the usual parton distribution sets, each 
parton distribution presents its own singularities. As an 
example, in the MRST2001 parametrization |^, we have 

xq{x, Qo) = ^(1 + + Cx){l — 


with Esea. = —0.26, = —0.33, = 0.09. In fact, 

these singularities do not correspond to any singularity 
present in hadronic cross sections and, conversely, cross 
section singularities are not present in parton distribu¬ 
tions. There must therefore exist a mechanism through 
which the singularities in partonic distributions disap¬ 
pear and cross section singularities arise when goes 
to zero. Such a mechanism is unknown and seems for¬ 
bidden by Regge theory. In this framework, a singular¬ 
ity structure common both to parton distibutions and to 
hadronic cross sections is the most natural choice. 

At that level, one may ask whether the Regge fit of 
Pi is compatible with pQCD and whether it is possible 
to have the same singularities in allparton distributions. 
Actually, although Regge theory PP and DGLAP P 
evolution both provide well-known descriptions of the 
structure functions P 0, 0, P, , the con¬ 

nection between the two approaches is unclear. In this 
paper, we will confront the triple-pole parametrisation. 
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for each parton distribution, and evolve it with DGLAP. 
This is done by fixing the initial distribution at Qq in or¬ 
der to reproduce the F^ value obtained from the global 
QGD fit p. We shall see that we are able to produce a 
fit to experimental data which is compatible both with 
Regge theory and with the DGLAP equation. This com¬ 
parison of two aspects of the theory will allow us to split 
the F 2 structure function in smaller contributions and 
to predict the density of gluons, which is generally not 
accessible directly from Regge fits. 


Varying the initial scale Qq, we can predict where per¬ 
turbative QGD breaks down. However, due to the appli¬ 
cation domain o of the global fit, the Regge constraint 
on the initial parton distributions is not valid at large x. 
In order to solve this problem, we use the GRV98 parton 
distributions 0 at large x {x > ccgrv)- We shall argue 
that the results do not significantly depend on the choice 
of the large x parametrisation. Since we will use leading 
order (LO) DGLAP evolution, one can choose any of the 
usual PDF sets to extend our fits to large x. 


We will show that, within a reasonable region of Qq 
and XGRV, the triple-pole pomeron model provides an 
initial condition for LO DGLAP evolution which repro¬ 
duces the experimental data. The scale Qg should be 
considered as the minimal scale where perturbative QGD 
can be applied. 


Since a good precision on the gluon density is of pri¬ 
mary importance for the LHG, it is also very interesting 
to look at the prediction of this model for the density of 
gluons. We will see that the densities we obtain are of 
the same order of magnitude as in the usual DGLAP fits. 


One should mention that such an extension of the triple 
pole Regge fit_^ a DGLAP evolution has already been 
introduced in [^. However, as we will see, our approach 
here is different: our parametrisation is much more con¬ 
strained, we are able to extract a gluon distribution and 
all the distributions have the same singularity structure. 
There are also some less important differences in the 
treatment of the large-a: domain. 
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II. PERTURBATIVE QCD AND REGGE 
THEORY 


A. Perturbative QCD 


In pQCD, the high behaviour of Deep Inelastic 
Scattering (DIS) is given by the DGLAP evolution equa¬ 
tions 1^. These equations introduce the parton distri¬ 
bution functions qi{x, Q^), qi{x, Q^) and g{x, Q^), which 
represent the probability of finding, in the proton, respec¬ 
tively a quark, an anti-quark or a gluon with virtuality 
less than and with longitudinal momentum fraction 
X. When oo, the evolution of these densities 

(at fixed x) are given by the DGLAP equations 


Q^dQ2 


/qi{x,Q^)\ 

qi{x,Q'^) 


27r A ^ 


I ■ At 93 

\ ^g<! Pgq 



( 2 ) 
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at leading order. Using these definitions, we have 


following parametrisation for the pomeron term 
a((3^)ln^ [ly/iyoiQ^)] +c{Q^), 
corresponding to a triple pole in j-plane Q 
a 2aln(r'o) aIn^(r'o)-I-c 


(4) 


(J - 1)^ (J - 1)^ 


J-1 


This seems to be the preferred phenomenological choice 
at = 0 Note that the upper expression, given in 
terms of ly and Q^, can be rewritten in terms of and 
X = Q^l{2iy). 

In a previous paper P|, we have also shown from 
unitarity constraints that we can extend the Gribov- 
Pomeranchuk argument about factorisation of residues 
to any number of thresholds and to any type of singular¬ 
ities. Hence, if we parametrise the pp and the 'y^*'ip cross 
sections, we can predict the cross-section using 

the t-channel unitarity (tCU) relation 


QlJ Q 2 ) 


^7p(j> Ql)^7p(j! Q2) 

^PP (J ) 


-|-finite terms 


F2{x,Q‘^) =xY^el. [qi{x,Q^) -f qi^x^Q"^)] , (3) 

i 

where the sum runs over all quark flavours. 

The usual way to use this equation is to choose a set 
of initial distributions 93 (x, Qq, a) to compute 93 (x, a) 
using @ and to adjust the parameters a in order to re¬ 
produce the experimental data. This approach have al¬ 
ready been successfully applied many times 0 , 0 , 0 , 0)01 
and is often considered a very good test of pQGD. Nev¬ 
ertheless, these studies do not care about the singularity 
structure of the initial distributions, ending up with re¬ 
sults that disagree with Regge theory, and presumably 
with QGD. 

B. Regge theory 

Beside the predictions of pQGD, we can study DIS 
through its analytical properties. In Regge theory 00 , 
we consider amplitudes A{j, t) in the complex angular 
momentum space by performing a Sommerfeld-Watson 
transform. In that formalism, we choose a singularity 
structure in the j-plane for the amplitudes. The residues 
of the singularities are functions of t and this technique 
can be applied to the domain cos(0t) ^ 1- For example, 
we can fit the DIS data or the photon structure function 
at large ly (small x), and the total cross sections at large 
s. 

The models based on Regge theory [H HI 0 0 
use a pomeron term, reproducing the rise of the struc¬ 
ture function (cross sections) at small x (at large s), and 
reggeon contributions coming from the exchange of me¬ 
son trajectories (a and /). We shall consider here the 


for the amplitudes in the j-plane. This relation proves 
the universality of the singularities, in other words, all 
singularities present in yp interactions also appear in 77 
interactions. We have successfully applied the tCU rules 
to the case of double and triple pole pomeron models in 
the region O- Therefore, since our fit keeps consistency 
with the fits in 0 for < Qq, we also be used to repro¬ 
duce the experimental results for Qf < Q 2 £ Qo- 


III. INITIAL DISTRIBUTIONS 

In our approach it is not possible to dissociate the soft 
sing ularity from the perturbative ones, as was done in 
[13 ■ However, it is possible to assume that the pertur¬ 
bative essential singularity (at j = 1) becomes a triple 
pole at small Q^. This may come from further resum¬ 
mation of pQCD 0. We thus have two regimes: for 
Q > Qq, we have a perturbative DGLAP evolution with 
an essential singularity, while for < Qq, the Regge fit 
applies, and F 2 behaves like a triple pole at small x. 

Due to the fact that the domain o does not extend 
up to X = I, we have used the GRV98 0 parametrisa¬ 
tion at large x, i.e. for x > xgrv- R is worth mention¬ 
ing that, in the DGLAP equation the evolution for 
X > Xgrv does not depend on the distributions below 
XGRV- This means that the evolution of the GRV98 dis¬ 
tribution functions for x > xgrv is not influenced by the 
parametrisation we will impose for x < xgrv- 

Since we want to use our fit 0 to F 2 , for X < XGRV, 
we want to have an initial distribution of the form {Qq is 
the scale at which we start the DGLAP evolution) 

F 2 {x, QI) = olog^(l/x) -I- 61og(I/x) -I- c-l- dx’', (5) 
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i.e. described by a triple pole pomeron and an /,a 2 - 
reggeon trajectory (77 = 0.31 as given in Q)- Once 
we have that initial distribution, we can evolve it with 
DGLAP and compare it with experimental data. 

However, the DGLAP equation m does not allow us to 
compute F 2 directly. Performing linear combinations in 
one can easily check that the minimal set of densities 
needed to obtain F 2 from the DGLAP equation is given 

by 


Moreover, we want to fix F 2 {Qq) to be equal to F^^'^ 
obtained from our previous global fit (each quantity with 
a superscript refers to the corresponding quantity ob¬ 
tained from the Regge fit in P|). Since T is entirely 
known, this constraint fixes all the S parameters through 
the relation 


<Pe = - z - 


(j) = a,b, c, d. 


(9) 


T = X [(u+ -b c+ -f t+) - (d+ -b s+ -b 6 +)] , ( 6 ) 
E = X [(u+ -b c+ -b t+) + {d+ -b s+ -b b+)] , (7) 
G = xg, ( 8 ) 


where q'^ = q + q ior q = u,d, s,c,t,b. The evolution 
equations for these distributions turn out to be 


(5^5q2T(x, = 


27r 


. e 




Q^dQ2 



CG f^^fPgq 2nfPqg\ 

Jx V £"? ^39 ) \^) 


and F 2 is then given by 


^ 5E + 3T 

ro, = -. 

18 

This clearly shows that, if we want to use as the ini¬ 
tial condition for a DGLAP evolution, we need to split 
it into T and S contributions, but we also need to intro¬ 
duce a gluon density. In this way, using 0 as the initial 
condition for the evolution allows us to predict the gluon 
distribution function. 

Since, below Qo; bo not use singularities of order 
larger than 3, we expect this behaviour to be valid for the 
T and S distributions. The natural way of separating the 
initial F 2 value given by 0 is thus to consider both T 
and E as a sum of a triple pole pomeron and a reggeon. 
The gluon distribution, being coupled to S, should not 
contain any singularities either. Thus, we can write 

r(x, Qo) = ot log^(l/a;)-b log(l/a;)-b ct- b dra;’’, 
S(x, Qo) = as log^(l/a;)-b 6 e log(l/a;)-b cs-b dsa:'', 
G{x,Ql) = aG^og^{l/x) + bG\og{l/x) + CG + dGX^. 

Most of the 12 parameters in these expressions are 
constrained. First of all, since the triple-pole pomeron, 
describing the high-energy interactions, has the vacuum 
quantum numbers, it will not be sensitive to the quark 
flavours. This means that, at high energy, one expects 
T ^ 0. Therefore, we set ar = bx = ct = 0. Then, since 
we connect our parametrisation with GRV’s at xqrVi we 
want the distribution functions to be continuous over the 
whole X range. Continuity of the T distribution fixes dx 
and we finally have 


Fix, Ql) = (xGRv, Ql) {" . 

V a:GRV / 


At this level, only the gluon distribution parameters 
are free. However, since the reggeon trajectory is ex¬ 
pected to be mainly constituted of quarks, we may ex¬ 
clude its contribution from the gluon density. Thus, we 
take dG = 0. Finally, we used continuity of the gluon 
density with the GRV distribution at xqrv to fix cg- 
We are finally left with only 2 free parameters: og and 
bG- 

Before going to the fits, one must stress that the GRV 
parametrisation at large x does not modify the triple pole 
singularity structure of the initial distributions. Actually, 
one may write the Mellin transform 

/*! /*a^GRV 

/ dx x^~^q{x) = / dx x^~^q^egge{x) 

Jo Jo 

+ / dxx^~^q^rv{x). 

J aJGRV 

In this expression, the hrst term generates the triple pole 
pomeron and the reggeon. The singularities of the second 
term come from the behaviour near x = 1. Since parton 
distributions behave like — cc)", the GRV parton 

distributions give the following contribution 

^ j + e + k-1- 

and the zeroes of the numerator cancel those from the 
denominator. Thus, using GRV at large x does not in¬ 
terfere with the singularity structure imposed from the 
low X parametrisation. 


IV. FIT 

We will fit the DIS data coming from Hl ll^ 1 ^ . 
ZEUSl2| m, BCDMSlH, E665[23, NMCjig and 
SLACj^. We will only consider data for F^. We have 
not included data from F'^^, Drell-Yan proccesses 

and Tevatron Jets for the following reasons 

• for many experiments, most of the data points are 
at large x or at low Q^. Thus, they do not constrain 
our ht much. 

• some experiments allow the determination of the 
valence quark distributions. We do not need them 
here since we only want the T, S and gluon distri¬ 
butions. 
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Since we want to test the domain common to Regge 
theory and to the DGLAP evolution, we only keep the 
experimental points verifying 

fcos(di) < 

< < g2 < 3000 GeV^ (10) 

I a; < XGRv- 


We have tried several values of the initial scale Qq around 
5 GeV^. Given an initial scale, the Regge limit on cos(0t) 
translates into a natural value for xqrv 


„(o) 

2 ^GRV - 49 


( 11 ) 



A graph of that limit is presented in Fig. ^ However, 



FIG. 1: Natural value of xgrv as a function of the scale 

as one can see from Fig. |2 if we take that limit on x, we 
cut most of the high experimental points which are 
at large x. It is therefore interesting to extrapolate the 
initial distributions to larger x, and we have tried some 
higher values for ccgrv- 

V. RESULTS 

The results of the fits are given in tabled as a function 
of Qq and ccgrv- We can see that this 2-parameters fit 
reproduces very well the experimental points in (IIOII for 
Qq > 3 GeV^ and xgbv <0.1. The values of the fitted 
parameters, as well as the constrained parameters are 
given in tableim 

We show the initial distributions and the Ff plot for 
Qq = 5 GeV^ and xgrv = 0.1 in Fig. |2|and Figs. I5I6I 
respectively. 

In Fig. 21 we have compared the gluon distribution ob¬ 
tained from our fit with that of the well known DGLAP 


FIG. 2: Experimental points, Regge domain limit and fit do- 
main limits for Qo = 5 GeV and a;GRV = or 0.1. It 

clearly appears that without extrapolation, we miss the high- 
points. 



X 


FIG. 3: Initial distributions for Qq = 5 GeV^ and xgrv = 0.1. 
'? 2/3 = x{u'*' -t c+ -I- t+) and g_i /3 = a;(d+ -t s+ -I- 6+) 

fits like GRV0, CTEQH and MRST0. One can see 
that our gluon distribution is of the same order of mag¬ 
nitude as that from other DGLAP fits. 

It is also interesting to check whether our results de¬ 
pend on the choice of the large x parametrisation. Since 
the DGLAP evolution equation couples the small x dis¬ 
tributions to the large x ones, at first sight, our results 
may depend on such a choice. However, looking at the 
studies of the PDF uncertainties, it can be seen that the 
large x behaviour of the T and E distributions hardly 
depends on the chosen fit down to x « 0.1. Moreover, in 
the large x limit, the splitting matrix can be written 

(P,, P.o] ^ 1 \ 

\Pm P 99 I ( 1 -^)+ V ■ ’ 
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FIG. 4: Fitted gluon distribution compared with some well 
known parton distributions 


XGRV 

‘^GRV 

0.1 

0.2 

Ql 


n 

X^/n 


n 

X^/n 

X^ 

n 

X^/n 

10.0 

484 

515 

0.939 

561 

581 

0.966 

774 

639 

1.212 

5.0 

557 

577 

0.966 

676 

686 

0.985 

862 

744 

1.159 

3.0 

633 

591 

1.071 

741 

735 

1.008 

- 

- 

- 


TABLE I: for various values of Qq and xgrv- (n is the 

number of experimental point satisfying nnn 


Thus, in the large x region, the gluon distribution and the 
sea are not coupled. Since, in our method, both T and E 
are fixed, we study the influence of the gluon distribution 
on F 2 ■ Due to the fact that these are not coupled at large 
x, we expect that our fit does not depend on the large x 
behaviour of the distributions. 


Qo 

3.0 

5.0 

10.0 

CL^p 

0.00541 

0.00644 

0.00780 

h^p 

0.0712 

0.0990 

0.142 

C~fp 

0.00541 

0.0064 

0.00780 

d'yp 

0.890 

1.06 

1.27 

xgrv 

‘^GRV 

0.1 

■^GRV 

0.1 

•^GRV 

0.1 

(It 

-0.0722 

0.167 

-0.0478 

0.166 

0.0101 

0.165 

ac 

0.147 

0.00617 

0.0908 

0.0271 

0.158 

0.131 

ha 

- 0.852 

0.718 

0.193 

0.822 

0.178 

0.419 

CG 

3.45 

-0.495 

0.595 

-0.851 

0.0299 

-0.463 


TABLE II; Values of the parameters for < Qq < 10 GeV^ 
and a^GRV < 0.1. Only aa and bo are fitted, while the other 
parameters are constrained. 

Furthermore, one can see that the of the fit remains 
of order 1 for 0.04 < xgrv ^ 0.15 and grows when we 
take a;GRV 0.01 or smaller. At that point , parton dis¬ 
tributions depend on the chosen parametrisation and the 
one we used, namely GRV98, does not take into account 
the latest HERA points. If we want to go to smaller val¬ 
ues of xgrVj we need a more recent parametrisation and 
thus a NLO study. Note that the interval on xgrv for 
which we have an acceptable hardly depends on Qo 
for Qq in [3,15] GeV^, and that the of the fit does not 
change very much in that domain. 

Unfortunately, it is quite hard to determine a unique 
scale Qq or xgrv from the fit. From Tabled it is clear 
that CCGRV can be taken to be 0.1 but can not be pushed 
up to 0.2. But, as we have argued, for such values of Qq 
and such high a^GRV, we are outside the domain o and 
we may not ensure that Regge theory will still be valid at 
X = 0.1 and Q^ = Qq- We can thus adopt two different 
points of view: 

1 . we stay in the domain O- We have thus xgrv = 

^GRV Qo down to 3 GeV^. The 

problem is that as Qq goes down, xgrv goes down 
too. And, since high Q^ experimental points have 
large x values, we do not test pQCD over a large 
range of Q^ values. This effect can be clearly seen 
in Fig. [3 where we have plotted the experimental 
points, the Regge domain limit and the fit domain 

for Qq = 5 GeV^ and xgrv = ^grv f'-f- 
therefore difficult to predict a “best value” for Qq. 

2. we extrapolate the Regge fit outside the do¬ 
main O- The amount of points concerned by this 
interpolation can be seen in Fig. [3 In such a 
case, depending on our confidence in this extrap¬ 
olation, we can consider that pQCD applies down 
to 3 GeV^ or 5 GeV^ and xgrv « 0.1. This value 
is compatible with the HERA predictions as well 
as with the Donnachie-Landshoff prediction [l^ . 
Below 2 GeV^, the is larger than 1 whatever 
3^GRV is and values of the initial scale smaller than 
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1.45 GeV^ are ruled out at the 90% confidence level. 


VI. CONCLUSION 

In this paper, we have shown that it is possible to use a 
very simple analytic form, namely a triple-pole pomeron 
and a reggeon, as an initial condition for the DGLAP 
evolution. Applying the constraint from a global QCD 
fit obtained in a previous paper P| as well as some ex¬ 
pected properties of the parton distribution functions, we 
have shown that we can fit the DIS data in the domain 
Qq > 3 GeV a: < 0.1 and cos(0t) > 49/(2mp). This 
fit has only 2 free parameters in the gluon distribution. 

Our fit is at the interplay between Regge theory and 
pQCD. We have thus proven that Regge theory can be 
used to extend a pQCD evolution down to the non- 
perturbative domain. From the fit, we can also say that 
the scale down to which we can apply pQCD is of the 
order of 3-5 GeV^. 

Moreover, we have seen that our approach can be used 
to split F 2 in T and E-components with precise physical 
properties. In this way, it is of prime importance to point 
out that all the initial distributions have the same sin¬ 


gularity structure, which is rarely the case for the usual 
parton sets. Since E is coupled to the gluon distribution, 
the latter can also be predicted. We have shown that 
the fitted gluon distribution is of the same order of mag¬ 
nitude as the gluon distributions obtained by the usual 
DGLAP fits to DIS data. 

By requiring the same singularities in each distribu¬ 
tion, we have seen that we are able to construct a full 
model both for DGLAP evolution and Regge theory in 
the case of a triple-pole pomeron model. It should be in¬ 
teresting, in the future, to test if we can apply the same 
method to the case of double pole pomeron or Donnachie- 
Landshoff two-pomeron model. 

In the future, it should also be interesting to see if it 
is possible to adapt this point of view, in order to derive 
the triple pole pomeron form factors at high Q^. 
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